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In many statistical multifragmentation models the volume available to the 
nonoverlapping fragments forming a given partition is a basic ingredient 
serving to the simplification of the density of states formula. One there- 
fore needs accurate techniques for calculating this quantity. While the direct 
. Monte-Carlo procedure consisting of randomly generating the fragments into 

f-*) \ the freeze-out volume and counting the events with no overlapped fragments 

CN \ is numerically affordable only for partitions with small N , the present paper 

' proposes a Metropolis - type simulation which allows accurate evaluations 

Q ■ of the free volume even for cases with large N . This procedure is used for 

. calculating the available volume for various situations. Though globally this 

T-H \ quantity has an exponential dependence on A^, variations of orders of mag- 

^ \ nitude for partitions with the same may be identified. A parametrization 

^ ' based on the virial approximation adjusted with a calibration function, de- 

■ scribing very well the variations of the free volume for different partitions 
. having the same N is proposed. This parametrization was successfully tested 
\ within the microcanonical multifragmentation model from [Al. H. Raduta 

3 : and Ad. R. Raduta, Phys. Rev. C 55, 1344 (1997); ihid., 56, 2059 (1997)]. 

' Finally, it is proven that parametrizations of the free volume solely depen- 

■ dent on are rather inadequate for multifragmentation studies producing 
I . important deviations from the exact results. 

^ ■ PACS number(s): 25.70.Pq, 24.10.Pa 

c : 

> : 

X ■ I- INTRODUCTION 

M': 

Nuclear multifragmentation is intensely studied both theoretically and experimentally 
from more than fifteen years. Since it concerns the decay of excited systems formed in 
violent heavy ion collisions supposed to be statistically equilibrated, this process has often 
been described by means of statistical models (e.g. f^-^). A basic ingredient of some of these 
models is the volume disponible to nonoverlaping fragments (often called free volume 
and denoted herein by ^eff), positioned into a freeze-out volume Q. This serves to the 
simplification of the density of states formula and implicitly to an easier evaluation of the 
various statistical observables. For the case of A^ noninteracting particles, ^eff is simply 
given by When the fragments are not allowed to overlap each other - which is the 
physical case, this quantity is drastically limited. In statistical models Qeff manifests as 
a factor entering the statistical weights of the fragment partitions and thus the accuracy 
of the model predictions is directly influenced by the quality of the ^eff parametrization. 
While an exact analytical evaluation of the A^-dimensional integral defining fleff is up to 
now practically impossible, various approximations for this quantity have been reported 
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since the beginning of the multifragmentation studies [§|,||,|T1. 

The topic of the effective volume goes far beyond the nuclear fragmentation studies being 
a fundamental thermodynamical problem. Fields like chemical physics, solid state physics 
and others are dealing with such problems. Here it is worth mentioning the approaches 
corresponding to infinite systems composed of identical or non identical spherical fragments 
proposed more then thirty years ago 0,0. While the above approaches are accurate for 
the systems they have been designed for, from the point of view of the nuclear multifrag- 
mentation studies they have two important drawbacks: they do not account for the finite 
volume of the freeze-out recipient in which fragmentation is supposed to take place and they 
have complicated expressions making them unsuitable to be employed in a time consuming 
Metropolis sharp microcanonical multifragmentation calculation. 

So far, accurate predictions of ^eff have been performed via a simple Monte Carlo 
method consisting of randomly generating the fragments into the freeze-out volume and 
counting the events in which no two fragments are overlapped However, though it 

provides exact evaluations of ^eff (within the statistical error), this method is numerically 
affordable only for cases with small number of fragments. Indeed, since ^eff is exponentially 
decreasing with A^, accurate evaluations for large are possible only with prohibitively large 
numbers of events. 

The present paper proposes a Metropolis-type simulation procedure which makes possible 
exact evaluations of ^eff even for partitions with large number of fragments. Since variations 
of X (defined as Qeff/^^) of orders of magnitude are observed for partitions with the same 

but different fragments size distribution a parametrization of x only dependent of N 
appears to be inadequate. For this reason a parametrization based on an adjusted virial 
approximation is proposed. 

The paper is organized as follows: Section |T| gives a description of the standard Monte 
Carlo method and proposes a Metropolis type procedure allowing exact calculations of x 
even for large N. Section presents a brief overall view on the results of the proposed 
method. In Section a virial parametrization of x is proposed. Applications and tests of 



the resulted parametrizations within the microcanonical multifragmentation model ||TO||ri 
are presented in Section |V|. The conclusions are drawn in Section |VT|. 



II. CALCULATING PROCEDURE 

Suppose that we are dealing with N spherical fragments of sizes Ai, i = 1, . . . , N resulted 
from a spherical source nucleus of size A: Y^iLi = ^) positioned at rj , i = 1, . . . , into 
a freeze-out recipient having the volume Q = n Qq, where Qq is the volume of the source 
nucleus at normal nuclear matter density: Qq = 47rrQy4/3, tq = 1.2 fm. The fragments are 
not allowed to overlap one another or the recipient "walls" meaning that they are subject 
to the constrains: Uf<j % n,=i with % = Qdvi - rj \ - {Ri + Rj)), 9^ = e{R - Ri - \ri\), 

where R = roA^/^n^/^, Ri = tqA^^ and G is the step function. The volume available to the 
A^ fragments writes: 

N N 

^e//(iV) = n/ dr, ^[n%- (2-1) 

1=1 ''^ i<j 
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From now on we will explicitely deal with the ratio between fleff and fl^ , x- As mentioned 
earlier, one can evaluate this quantity by means of a simple Monte Carlo technique consisting 
of placing randomly each of the N fragments in the considered freeze-out recipient. The 
event is considered successful if no two fragments are overlapping each other and no fragment 
is intersecting the recipient walls. Supposing that one obtains Ns successful events from Nt 
attempts, one can estimate x according to: 

X=^- (2.2) 



AT, 



Since the relative statistical error in evaluating is of the order of Xj^fW^^ it results 
that for accurate evaluations of small x values (which is the case of large A^) one needs a 
very large number of generated events [Nt) which makes the method unpracticable. For this 
reason, a Metropolis-type simulating procedure for calculating x is further proposed. 
Using (|2.1|) one obtains 

'''^^^''^ -^mh. (2.3) 



and, by recurrence. 



N 

X=Vi9k fk, (2.4) 

k=l 

where gk = dr^ 91/Q = [{R - Rk)/Rf and 

r ^ In dri e[ dT2 ei ^12 • • • dr, 9,^ ■ ■ ■ 9k-i k , . 

dn 91 dr2 9^, 9u... In dn-i 9l_, 9, . . . 9,_, dr^ ^ ' ^ 

The last expression can be regarded from the point of view of a statistical ensemble cor- 
responding to a system consisting of the first k (g {2, . . . , A^} - we don't need to consider 
k = 1 since /i = 1) fragments from the considered partition, composed of configurations 
C : {rj, i = 1, k} and subject to the constrains nf</ % Ili=i This way, equation ( p.5| ) 
can be re-written as follows: 

/.= ^^^^l--f-^^ = (^i....^.-i (2.6) 

where (■)^ has the meaning of average value over the states of the above mentioned ensemble. 
Therefore, it is sufficient that one knows the values of fk for A; = 2, . . . , in order to obtain 
X- These values can be evaluated using ( p.6| ) within a Metropolis-type simulation which is 
further described. 

For evaluating one just has to generate a Markovian walk in the configuration space 
of the k fragment system according to the detailed balance principle: 

AC WcP{C C) = AC Wc'P{C ^ C), (2.7) 

where AC and AC are the elementary volumes of the configurations C and C, Wq and 
Wc' their statistical weights (here equal to unity) and P{C C) and P{C ^ C) are the 
probabilities of passing from C to C and respectively of making the reverse move. 
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The random walk is generated as follows. Suppose that the current state of the system 
is C. One chooses at random one of the k fragments, indexed by i and positioned at r,. 
One randomly re-positionates this fragment into a spherical volume having the same center 
as the freeze-out recipient and the radius R — Ri (this ensures the non-overlapping between 
the fragment and the recipient wall). We denote the chosen position of the fragment i with 
r^. At this point one has to check whether the constrains of the system are satisfied. So, 
if i k and the fragment in its new position is overlapping at least one of the remaining 
fragments indexed from 1 to — 1, the move is aborted and configuration C is reconsidered. 
Otherwise, the configuration C is correct and is considered as a new configuration of the 
system. The probability of this move is: 

The probability of the reverse move, similarly generated, writes: 

Taking into account that AC = drjllj^idrj, AC = dr-Ilj^idrj and that Wc = Wc, 
one may easily see that the detailed balance equation, ( p.7|) is satisfied. At this point the 
exploration of the configuration space is completely described. For calculating one just 
has to apply formula (|2.6|) , i.e. to determine the average value of the observable Ilfj/ (^ik over 
the states selected by the simulation. For a number of Ng successful (for which Yl^zl Oik = 1) 
accepted events the relative statistical error in estimating fk is of the order of l/^/Wl. 
Provided that fk was evaluated for each G {2, . . . , A^} using the above described sim- 



ulation, X can be obtained using (2]^). Assuming that each simulation (corresponding to a 
given k) is performed using successful accepted events, then the relative statistical error 
in evaluating x is of the order of -y/iV/iVg. 



III. OVERALL RESULTS 

In order to verify the accuracy of the proposed simulation we present in Fig. |1| a com- 
parison between x{^) calculated with the present simulation and that calculated with the 
direct Monte Carlo procedure briefly described at the beginning of Section I in the case of 
A = 100, f2 = 10 fio (i-6. n = 10). For this evaluation we considered A^ in a relatively 
small range (2 to 17) because for larger A^ the direct Monte Carlo calculation would require 
unreasonable computing time. The fragment partitions corresponding to each A^ are ran- 
domlyQ selected. This calculation was performed for A^s = 10^ successful events implying 



Starting from the source nucleus one randomly splits it in two fragments. Then, one randomly 
choses one of the resulted fragments and splits it randomly in two fragments. The process is 
repeated until fragments are obtained. This method has the advantage of generating partitions 
with X taking values from a wide range for a given N. The dispersions of x obtained using this 
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that the relative errors are smaller than the dimension of the points. The same Ng is used 
in all the calculations presented in this paper. One can observe that the results of the two 
calculations are practically identical. 

For having a preliminary overall look on the simulation results, we calculate x for ^ = 200 
and n =4, 6, 8, 10, 14, 20 for ranging from 2 to 30 with a step of 1. For each considered 
A^, 5 simulations are performed corresponding to randomly chosen partitions. This allows 
to evaluate the variation of x for partitions having the same N. The results are presented in 
Fig. ^. One can observe that while globally log^g X has a linear behavior, variations of orders 
of magnitude, strongly increasing with decreasing n can be observed for % corresponding to 
different partitions with the same A^. 

In order to test to what extent the assumption of nonoverlapping between the fragments 
and the recipients' wall is contributing to the above mentioned result, we run the same 
simulation without imposing this boundary any longer (i.e. the constrains are Ylf^jOij). 
Calculations of x a^re performed for A = 200 and n =3, 4, 6, 8, 10, 14, 20 for ranging from 
2 to 30 with a step of 1 and 10 randomly chosen partitions per each A^ value (see Due to 
this relaxation, the effective free volume is now larger so the dependencies logj^g x(^) are less 
abrupt than in the previous case (see Fig. |^) and their dispersion is smaller. Nevertheless, 
variations of x of orders of magnitude corresponding to partition with the same A^, strongly 
increasing with decreasing n, can be identified here as well. 

This suggests that parametrizations of x solely dependent on A^ are quite inadequate 
for cases with relatively small values of n - typically used in multifragmentation studies. 
One therefore needs a parametrization of x valid over a large range of A^, independent of A 
and taking into account the above-evidenced strong variation of x when A^ and n are kept 
constant. This problem is addressed in the next section. 

From now on we will concentrate exclusively on the initial boundary assumption (i.e. 
nonoverlapping between the fragments and the recipient walls). 



IV. PARAMETRIZATION 

The aim of this section is to provide an accurate parametrization of x which takes 
into account the strong variation of this quantity for partitions with the same number of 
fragments. Replacing the expression of x (corresponding to the case in which the fragments 
are allowed to intersect the recipients' walls): 

(4.1) 

r=i In^^i 

by the following "factorization" : 



generating procedure are even larger than those obtained using partitions realistically generated 
(i.e. by means of a microcanonical multifragmentation model). This can be easily checked by 
comparing the dispersions of x (corresponding to n = 6) from Fig. 2 with those from Fig. 8. 
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where 9ij = 1 — % and Pjj = {J^dvi J^drj dr^ /j^ dr^), one easily recognizes the 

two-body approximation due to Cole, Heuer and Charvet [§]. Note that here Pij has the 
meaning of probability of overlapping between the particles i and j when they are randomly 
generated into the volume Q. A more convenient form for (|4.2| ) is: 

N 

lnx' = T.^n{l-P,,) (4.3) 

i<j 

which for small values of Pij can be approximated by: 

N 

lnx. = -E^ir (4-4) 

i<j 

The last expression is the so called virial approximation proposed by Randrup, Robinson 
and Sneppen 0. While the exact expression of the two particle overlapping probability, Pij, 
was deduced in we prefer the simplest approach f^: 

Obviously, since it neglects the higher order interactions between fragments and allows in- 
tersections of the fragments with the recipients' walls, this approximation overestimates 
X- Nevertheless, apart from other approaches (see e.g. [^,|^ used in chemical physics cal- 
culations) the "virial" approximation has the advantage that it is very simple and that it 
explicitely accounts for the volume of the freeze-out recipient (thus being appropriate for 
multifragmentation studies). As we shall further see, this approach provides the basis for a 
very accurate parametrization of the free volume. 

In order to estimate the difference between the virial evaluations (xv) and the exact X) 
evaluated by means of the method proposed in Section we represented in Fig. ^ the ratio 
liO-X/ InXi; calculated over a large range of N (sufficient for the systems currently considered 
in multifragmentation studies), for four values of n: 4, 6, 8 and 10. For each considered N 
a fragment partition was generated as in Section |TT1[ The considered source is A = 300. 
Though the fragment partitions are generated as to induce an important dispersion in x (see 
Fig. 1^), Inx/ InXv has a smooth behavior. In order to have an estimate of the fluctuations 
of the Inx/ ^^Xv{N) dependency, we evaluate this quantity for n = 6, A = 300, ranging 
from 2 to 24 and 10 randomly chosen partitions for each value of A^. (see Fig. |^). Though 
the fluctuations are increasing with decreasing A^ they lay in reasonable limits and are 
practically negligible for A^ > 10. 

Therefore, one can obtain accurate fits of lnx/lnXt,(A^) with appropriate functions. As 
shown in Fig. |^, the function: 
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provides a good fit for In xl In Xv versus A^, for all considered n. The corresponding pa- 
rameters are listed in Table |. Knowing /(A^) for a given n, the real Inx can be expressed 
as: 

lnx = lnx. /(iV), (4.7) 



with Inx^ given by (4^). In order to test this parametrization, in Fig. ^ are represented 



versus x calculated using the method from Section II and x evaluated by means of ( [4.7| ) 
for n =4, 6, 8, 10. The fragment partitions were considered as in the previous paragraph. 
The range of iV is considered only up to = 100 just for the clarity of the plot. As one can 
see, the agreement is remarkably good. It thus appears that the present parametrization has 
the ability of describing accurately the strong variations of x fo^' different partitions having 
the same A^. Moreover, since for a given partition x only depends on the fragment relative 
size distribution, it results that this parametrization is independent on A. Some applications 
and tests of this parametrization are presented in the next section. 



V. APPLICATION TO A MICROCANONICAL MULTIFRAGMENTATION 

MODEL 

In this section some supplemental tests for the accuracy of the parametrization proposed 
in Section |I^ performed with the sharp microcanonical multifragmentation model from 
Ref. [|l37|nil are presented. This model sharply conserves the number of nucleons (A) the 



number of protons (Z), the total energy {E) and the total momentum (P) and assumes 
equal probability of appearance of all possible configurations C\ : {A,, Zj, r^, e^, pi, i = 
1, . . . , A^}, where the parameters under the brackets are respectively the mass number, the 
atomic number, the position, the excitation energy and the momentum of each of the A^ 
fragments composing the respective configuration. Integrating the expression of the total 
number of states of the system over the fragment momenta one reduces the configuration 
space and works into a new one, composed of configurations C : {Ai, Zi, r^, ei, i = 
1, . . . , A^} having statistical weights of the form: 



Wc = — TT 



{mAiY 



2n {2tiK) 



\N- 



r(|(Ar-l)) ^rnA)i 



(5.1) 



where pi is the level density of the fragment i, m is the mass of a nucleon, r() is the 
hyper-geometric function and K is the kinetic energy of the given partition and can be 
expressed as: K = E + J2i Bi — — HiKj^ij- Here Bi is the binding energy of the 
fragment i and Vij is the two-body Coulomb interaction energy. The weight of a configuration 
C : {Ai, Zi, e,, z = 1, . . . , A^} corresponding to a new space where the fragment positions 
are now missing, can be written as: 



N , N 



Wc' = t{l dr^PFc^rn^, (5.2) 

7_i -/ri ■ 



i=i-"-' i<j 

where fl, 9\ and 9ij have the same meaning as in the previous sections. Now, it can be 
observed that if one replaces the multi-body Coulomb interaction depending on the fragment 
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positions (entering in K and therefore in Wc) by a single particle mean-field approach, Wc 
nolonger depends on and (|5.2|) can be rewritten as: 



Wc = X Wc. (5.3) 

As one can see, in this case x nianifests as a factor entering the statistical weight of a 
given configuration, C. The average value of any system observable X is evaluable in this 
new ensemble by: {X) = J^c Wc'Xc / Y.c Wc- For doing this one applies identically the 
Metropolis simulation proposed in [|lO|,|lT| except that the positions are no longer gener- 



ated. The resulted correction factor, a is the same as in Ref. |IO| , p!T| except for the factor 
x(C^+i)/x(C'iv) appearing now due to the modification of the configuration weights (see 
equation. (|5.3|)). For the Coulomb interaction the Wigner-Seitz approach is employed ad- 
justed with a factor as to provide accurate descriptions of the total Coulomb energy even in 
cases with small A^: 

Vc = Vws g{N), (5.4) 

with 

where = Ri m}^^ , Vws is the Wigner-Seitz energy and the factor g{N) = (1.07675 + 
18.5266) / (A^+ 14.9439) correspond to = 6, which is the case for which the microcanonical 
calculations are here performed. This factor was evaluated by fitting the ratio between the 
average Coulomb energy (evaluated by uniformly generating the fragments from a given 
partition into the spherical recipient such as they are not overlapping each other - this 
is performed by simple Metropolis moves) and the corresponding Vws versus A^. It is 
worth mentioning that, included in the microcanonical model, this parametrization insures 
practically identical results with those obtained with the unmodified version of the model. 

For testing the accuracy of the free-volume parametrization two types of simulations 
are performed via the microcanonical multifragmentation model. In the first simulation we 
just replace the multi-body Coulomb interaction energy with the Wigner-Seitz approach ex- 
posed above. Except this detail, the simulation is identical with that from Refs. [rU,|rT| and 



therefore the fragment positions are explicitely generated at each step, the forbidden con- 
figurations (in which fragments are overlapping) being automatically rejected. The second 
simulation is that formulated in the previous paragraph and includes both the Wigner-Seitz 
parametrization and the free-volume parametrization. In principle, a correct parametriza- 
tion of X should provide closely identical results for any system observable. Here we choose 
to compare the mass distributions calculated within the two simulations for the source nu- 
cleus (70, 32) and the excitation energies -Eex=2 MeV/nucleon and -E'ex=6 MeV/nucleon. 
These are represented in Fig. ^. As one can see, the agreement is very good for both cases 
fact which supports the present parametrization. 

An interesting question is to what extent a free-volume parametrization only dependent 
on A^ is appropriate for nuclear multifragmentation studies. For answering this question 
one needs to obtain for each of the above considered parametrization of x only 

dependent on A^. To this aim, within the first of the two simulations mentioned in the 
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previous paragraph, for each situation [(70, 32), Eex='2, 6 MeV/nucleon] we select at random 
different fragment partitions, corresponding to each N appearing in the simulation, for 
which we calculate x- These are represented in Fig. |^. Since the obtained dependence is 
rather linear in logarithmic scale we fit these distributions by means of linear functions (see 
Fig. 1^. We then employ the obtained parametrizations x(A^) in the multifragmentation 
model and we perform the second type of simulation described in the previous paragraph. 
The results are represented in Fig. ^ by dashed-dotted lines. As one can observe, the 
difference between the mass distributions obtained with this last simulation and the mass 
distribution corresponding to the first simulation (mentioned in the previous paragraph) are 
quite significant. This confirms the fact that the dispersion in x for different partitions having 
the same N which was illustrated in Figs. ||, || produces indeed important deviations of the 
model results from the real values and, therefore, a parametrization of x solely dependent 
on N appears to be inadequate. This result is important since it is well known that most of 
the existent multifragmentation models employ parametrizations of x only dependent of N. 

Finally, a comparison between x{^) corresponding to the partitions considered in the 
previous paragraph (for the source (70, 32) and Eex=2, 6 MeV/nucleon) and the SMM 
parametrization is performed. According to the SMM prescription one has: x = {[1 + 
d{N^^^ — 1) / (rQA^^^)]^ — l}^/n^ where d=1.4 fm. The x{^) dependencies calculated as 
described in the previous paragraph are represented in Fig. |^ together with the SMM one 
calculated for n = 3 - the usual SMM assumption (since in the SMM model the free volume 
parametrization has no explicit dependence on n we prefer to use its standard freeze-out 
assumption, n = 3). Analyzing Fig. ^ one can observe that the SMM n = 3 curve is flatter 
than our calculations performed for the n = 6 case. Taking into account that the lnx{N) 
slope is increasing with decreasing n (see Fig. it follows that the differences between 
the SMM curve and the exact results corresponding to n = 3 are even larger than those 
observed in Fig. ^ between the exact calculations with n = 6 and the SMM n = 3 curve. 



VI. CONCLUSIONS 

In summary, a Metropolis-type procedure was designed in order to describe the volume 
available to N spherical fragments positioned into a spherical recipient. Apart from the 
direct Monte Carlo technique, this method allows one (for the first time) to make x evalua- 
tions even for large number of fragments. Calculations of x for randomly chosen partitions 
corresponding to various n and values are evidencing that while, for a given n, x has 
a more or less global linear dependence of A^ in logarithmic scale, variations of orders of 
magnitude may be identified for x calculated for different partitions corresponding to the 
same n and A^. This suggests that parametrizations of x only dependent on A^ may be in- 
adequate for multifragmentation studies. For this reason a parametrization based on virial 
approximation adjusted with a calibration function as to fit the exact x values is pro- 
posed further-on. The obtained parametrization describes the x variations corresponding to 
different partitions having the same A^ remarkably well. Calibration parameters have been 
evaluated for four values of the volume of the recipient, corresponding to n = 4, 6, 8, 10. 
The parametrization is further tested with the microcanonical multifragmentation model 
from Ref . • For doing this, the simulation is run in two forms: first with the inclusion 



of the Wigner-Seitz approach for the Coulomb interaction but maintaining the fragment 



9 



positioning into the freeze-out volume and the hard core repulsion and then with including 
both the Wigner-Seitz approach and the free-volume parametrization obtained in Section 
P^ . Mass distributions obtained with these two types of simulation are compared in two 
input situations. A very good agreement is to be noticed which supports once again the 
parametrization obtained in Section Finally a study is made concerning the extent to 
which a parametrization of x only dependent on is appropriate for multifragemtation 
studies. To this aim, the second simulation was run again this time using a parametrization 
of X only dependent on N . In this last case, important deviations between the mass yields 
obtained with the first simulation and those obtained with the second one are to be noticed. 
This suggests parametrizations of x only dependent on may induce important devia- 
tions of the models' results from the real ones. Comparisons between the evaluated x(A^) 
dependencies and the corresponding SMM parametrization show important discrepancies. 

Presently, there is a strong debate concerning the role of the freeze-out volume in de- 
termining the behavior of various thermodynamical quantities (such as the nuclear caloric 
curve). It was shown that the inclusion of a free volume parametrization in a statistical 
multifragmentation model is crucial for a good theoretical description of various experimen- 



tal isotopic caloric curves [12|. The value of the freeze-out volume was proven to influence 



the shape of the caloric curve and even the order of a possible phase transition: while larger 



volumes are generating caloric curves with continuously increasing plateau-like regions |11 



(reflected in positive peaks in the heat capacity curves), smaller freeze-out volumes seem to 
encourage the occurrence of back-bendings in the caloric curves, reflected in negative regions 
of the heat capacity curves |jl3|]. Even the definition of the freeze-out volume is still intensely 
discussed: should the statistical ensemble be isochore (see e.g. 0) or should it be isobar 
(which implies a smooth variation of the freeze-out volume with the excitation energy - as 
shown in Ref. this option leads to the appearance of back-bendings in the caloric curves) 
- see 0. Recently, nonspherical shapes of the freeze-out recipient have been investigated 
]r5| as well. Recent experimental evaluations of the freeze-out volume using two particles 
velocity correlations predicted values of n in a wide range: 2.5 - 12.5 [0. 

The author thanks Ad. R. Raduta for a critical read of the manuscript. 
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FIG. 1. Comparison between x values calculated by means of the direct Monte Carlo procedure 
and those calculated by means of the Metropolis-type simulation for A = 100, n = 10 and N 
ranging from 2 to 17. For each N a randomly generated partition was used. 
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FIG. 2. logioX versus iV calculated for six values of n and a source with ^ = 200. For each 
particular n and N 5 values of log^^Q x are represented, corresponding to randomly chosen partitions. 
The points corresponding to n = 6, 8, 10, 14, 20 are shifted upwards with the quantities specified 
in the legend in order to avoid overloading. 
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FIG. 3. Same as in Fig. |2| except that the calculations are performed without imposing the 
constrain of nonintersection between fragments and the recipient wall. Additionally, a calculation 
corresponding to n = 3 is included. For each particular n and N 10 randomly chosen partitions 
are used in the calculation. Points corresponding to various n are shifted upwards or downwards 
as shown in the legend in order to avoid overloading. 
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FIG. 4. Ratio between Inx and In^t; (see the text) versus the number of fragments (A^) calcu- 
lated for A = 300 and n = 4, 6, 8, 10 fitted with f{N) (see equation (|4.6|)). The corresponding 
fitting parameters are listed in Table |. 
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FIG. 5. lnx/^T^Xv{^) calculated for A = 300, n = 6 and 10 randomly chosen fragment partitions 
for each value of A^. 
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FIG. 6. Calculated lo 

§10 X versus corresponding to A = 300 and n = 4, 6, 8, 10. For each 
N, X was calculated for a randomly (see Section III) chosen partition. The circles correspond to 
calculations using the exact method from Section || and the triangles correspond to the adjusted 
virial approximation [eq. ( [4.7D ] using the parameters from Table |. 
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FIG. 7. Fragment mass multiplicities corresponding to the source (70,32) and the excitation en- 
ergies: 2.5 and 6 MeV/nucleon calculated using the microcanonical multifragmentation model. The 
continuous lines are calculated using the adjusted Wigner-Seitz (W.S.) approach (see the text) but 
keeping the fragment positioning into the freeze-out volume, the dashed lines are calculated using 



both the W.S. approach and the free volume (f.v.) parametrization (47) and the dashed-dotted 
lines correspond to the W.S. approach and a f.v. parametrization only dependent on N (see Fig. 
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FIG. 8. Inx versus AT calculated for the fragment partitions appearing in the fragmentation of 
the (70,32) source nucleus for the excitation energies 2.5 and 6 MeV/nuclcon. Linear fits of Inx(iV) 
are represented by dashed lines. The partitions corresponding to each N arc randomly selected 
from the output of the microcanonical model in which the W.S. approach was used. 
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FIG. 9. Inx(A^) calculated for the fragment partitions appearing in the fragmentation of the 
(70,32) source nucleus for the excitation energies 2.5 and 6 MeV/nucleon and Inx(iV) corresponding 
to the SMM free volume parametrization for the case of A = 70 and n = 3. 



n 

~ 

6 
8 

10 



2.189191 
3.665597 
0.927222 
1.117942 



1.096606 
-0.599464 
-1.519126 
-1.976695 



1.042162 
1.125311 
0.999927 
1.020913 



d 

-0.934725 
-1.202542 
-1.085739 
-1.100722 



0.711251 
0.432149 
0.863371 
0.870260 



TABLE I. The parameters of the function f{N) [([4. 61)] corresponding to four values of n 
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